Natural NMSSM Higgs Bosons 



S. F. King 1 ; M. Miihlleitner 2 1 R. Nevzorov 3 * K. Walz 2 § 

<^ 1 School of Physics and Astronomy, University of Southampton, 

(N Southampton, SO 17 IB J, U.K. 

> 

2 Institute for Theoretical Physics, Karlsruhe Institute of Technology, 
76128 Karlsruhe, Germany. 

^Institute for Theoretical and Experimental Physics, Moscow, 117218, Russia. 

Oh 

Abstract 

We study the phenomenology of Higgs bosons close to 126 GeV within the scale invariant 
unconstrained next-to- minimal supersymmetric Standard Model (NMSSM), focusing on the 
regions of parameter space favoured by low fine-tuning considerations, namely stop masses 
of order 400 GeV to 1 TcV and an effective \x parameter between 100-200 GeV, with large 
(but perturbative) A and low tan [3 =2-4. We perform scans over the above parameter space, 
[- — focusing on the observable Higgs cross sections into 77, WW, ZZ, bb, tt final states, and 

study the correlations between these observables. We show that the 77 signal strength may 
be enhanced up to a factor of about two not only due to the effect of singlet-doublet mixing, 
which occurs more often when the 126 GeV Higgs boson is the next-to- lightest CP-even 
t-H one, but also due to light stops (and to a lesser extent light chargino and charged Higgs 

loops). There may be also smaller enhancements in the Higgs decay channels into WW, ZZ, 
correlated with the 77 enhancement. However there is no such correlation observed involving 
the Higgs decay channels into bb, tt. The requirement of having perturbative couplings up 
to the GUT scale favours the interpretation of the 126 GeV Higgs boson as being the second 
lightest NMSSM CP-even state, which can decay into pairs of lighter neutralinos, CP-even or 
CP-odd Higgs bosons, leading to characteristic signatures of the NMSSM. In a non-negligible 
part of the parameter range the increase in the 77 rate is due to the superposition of rates 
from nearly degenerate Higgs bosons. Resolving these Higgs bosons would rule out the 
Standard Model, and provide evidence for the NMSSM. 
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1 Introduction 



The recent discovery of a new particle with a mass around ~ 125 GeV [l||2j is consistent with 
the Standard Model (SM) Higgs boson. In particular the observed decays and signal strengths 
into 77, WW, ZZ favour the interpretation that the particle is a neutral boson with spin-0. 
However more data is needed to assess its nature, and if careful studies of the signal strengths in 
different channels reveal discrepancies from the predictions of the SM then this would provide 
a window into new physics Beyond the Standard Model (BSM) |3j|. Supersymmetric (SUSY) 
models are a leading candidate for BSM physics and generically predict one or more light Higgs 
bosons whose properties may differ in detail from that of the SM Higgs boson. For example, if 
the cross section of Higgs production and decay into 77 were observed to be significantly higher 
than the SM Higgs prediction, then this could be due to the effects of SUSY particles in the 
loops [4-13 or suppressed couplings to 6-quarks leading to smaller total widths (9-24 . 



In the Minimal Supersymmetric Standard Model (MSSM) there exists an upper limit on the 
lightest Higgs boson mass of about 130-135 GeV, depending on the values of the parameters 



in the stop sector (see e.g. 25 and references therein). The MSSM can be consistent with a 
126 GeV SM-like Higgs boson in the decoupling limit. In this limit the lightest Higgs boson 
mass is given by 

m 2 h Rs Mf cos 2 2/3 + Am 2 h , (1.1) 

with the correction Am 2 being dominated by loops of heavy top quarks and top squarks. The 
ratio of the vacuum expectation values (VEVs) of the two Higgs doublets introduced in the 
MSSM Higgs sector is denoted by tan j3. In order to raise the Higgs boson mass to 126 GeV, 
we hence need at large values of tan /? a loop contribution of Am^ ~ 85 GeV which is nearly as 



large as the tree- level mass value. This leads to some degree of fine-tuning 26 



It has been known for some time that the fine-tuning of the MSSM could be ameliorated in 



the scale invariant Next-to-Minimal Supersymmetric Standard Model (NMSSM) [27,28 . With 



a 126 GeV Higgs boson, due to the fine-tuning of the MSSM, the NMSSM has emerged as a 



more natural alternative. In the NMSSM 29 -32] (for reviews see |33[|34] ) one singlet superfield 
S is added to the spectrum of the MSSM. The supersymmetric Higgs mass parameter fi is 
then generated dynamically through the coupling term XSHdH u . The upper mass bound of the 
lightest Higgs boson in the NMSSM becomes, 

\2 2 

m 2 h ps AfJ cos 2 2/3 + — - sin 2 2/3 + Am 2 h , (1.2) 

where v = 246 GeV. Contrary to the MSSM, for Xv > Mz, the tree-level contributions to m^ 
are maximized for moderate values of tan/3. For example, setting A = 0.6 and tan/3 = 2, these 
tree-level contributions raise the Higgs boson mass to about 100 GeV requiring Am/, ~ 75 GeV 
in order to match the 126 GeV Higgs mass value. The difference to the correction needed in the 
MSSM (numerically about 10 GeV) is significant as Am^ raises logarithmically with the stop 
masses and receives an important contribution from the stop mixing. 

In the NMSSM, depending on tan/3, A ~ 0.7 is the largest value in order not to spoil the 
validity of perturbation theory up to the GUT scale. The presence of additional extra matter, 
however, allows larger values of A to be achieved [35] . For example, adding three families of 
5 + 5 extra matter at a mass scale of 1 TeV increases the largest value to A ~ 0.8 for the same 
parameters as before. The above discussion shows that there is an argument from fine-tuning 
for extending the NMSSM to include extra matter. Such an NMSSM+ model with extra matter 
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has recently been discussed in [36]. 

In this paper we study the phenomenology of Higgs bosons in the mass range 124-127 GeV 
within the scale invariant NMSSM. To distinguish our study from the many NMSSM studies in 
the literature of a near 126 GeV Higgs boson, we shall focus exclusively on the regions of pa- 
rameter space favoured by low fine-tuning considerations, namely stop masses of order 400 GeV 
to 1 TeV and an effective /i parameter between 100-200 GeV, with large (but perturbative) A 
and low tan/3 =2-4. We shall allow for the possibility of extending the NMSSM to include 
extra matter as in the NMSSM+ |36|, so that A can be increased up to 0.8 at low energy scales, 
while remaining perturbative up to the GUT scale. We perform scans over the above parameter 
space, focusing on the observable Higgs cross sections into 77, WW, ZZ, bb, tt final states, 
and study the correlations between these observables. We show that the 77 signal strength 
may be enhanced up to a factor of two due to an enhancement of the production and/or the 
decay mechanism. While small stop mixing for light stops enhances the dominant production 
process through gluon fusion and suppresses the loop- mediated decay into photons, the latter 
can also be enhanced to some extent through light chargino and charged Higgs boson loops. 
Furthermore, the suppression of the dominant decay into bottom quarks due to singlet-doublet 
mixing entails a suppressed total width and hence enhances the branching ratio into photons. 
Since also the branching ratios into WW and ZZ are affected by such a suppression there is a 
strong correlation between these channels and the 77 enhancement. However there is no such 
correlation observed involving the Higgs decay channels into bb, tt which may be suppressed 
in the enhanced 77 region. The results include the possible presence of a second Higgs boson 
in the region of 126 GeV. The superposition of the rates of two nearly degenerate Higgs bosons 
also increases the event rate in the photon final state. Provided that such a signal can be 
disentangled from a singly produced Higgs boson in future, this will be a further strong test 
of the NMSSM. Our scan reveals that with the chosen small stop mass values it is difficult to 
get a lightest CP-even Higgs boson with mass around 126 GeV if in addition perturbativity 
constraints are imposed. It can only be achieved for large mixing values and inclusion of extra 
matter at ~ 1 TeV, while this is not necessary if the second lightest CP-even Higgs boson is 
demanded to have the same mass as the recently discovered new resonance. Interestingly, in 
this case the mass spectrum can be such that the heavier CP-even Higgs boson can decay into 
a pair of lightest neutralinos, CP-odd or CP-even Higgs bosons, which leads to distinctive final 
state signatures to be tested at the LHC. 

The work in this paper complements and goes beyond the other studies of a Higgs boson 
near 126 GeV in the NMSSM [H}jl^[l^[2l^|2^[2^[37Ht3] . For example the original observation 
that the di-photon channel may be enhanced due to strong singlet-doublet mixing due to the 
reduction of the bb partial width with a second lighter CP-even Higgs boson was made in [14] . 
This was followed by our proposal [19] of a set of benchmark points in which we studied, in 
addition to 77 and bb, also the channels WW and ZZ both for the case where the second CP- 
even Higgs boson is lighter or heavier than the 126 GeV one, for the case of light top squarks and 
gluinos. The channels were also studied in the framework of various versions of the constrained 

and a comparative study between the MSSM and 



NMSSM with relatively heavy stops 37 
NMSSM has been performed in [21] . In 



22 the effect of astrophysical and Dark Matter (DM) 



constraints on the lightest supersymmetric particle (LSP) in the NMSSM with a 126 GeV Higgs 
boson was taken into account with the main focus on the LSP being a singlino-like neutralino. 
Similar constraints were also applied to the constrained NMSSM with a dominantly Higgsino-like 



LSP 23 . The case of the Higgs boson mass spectrum in the complex NMSSM was considered 
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Ill 



38 leading to significant effects on Higgs phenomenology. In 39 scenarios were investigated 



where the two lightest NMSSM Higgs bosons are closely spaced near 126 GeV, leading to very 
enhanced decay rates, and in [24] scenarios with Higgs bosons both consistent with the LEP 
98 GeV excess and the 126 GeV Higgs boson of the LHC search were studied, while [l2j discussed 
the case where the NMSSM Higgs sector could both explain the 126 GeV discovery and the small 
excess observed by CMS at 136 GeV. In [40] a more complicated non-NMSSM model with extra 
singlets was proposed, while in [41] another alternative to the NMSSM involving singlet mass 
terms was studied. The case of fine-tuning in the NMSSM was analysed in [27] . Finally, the 
effects of combining the NMSSM with an inverse see-saw mechanism were considered in [42] . 

The layout of the remainder of this paper is as follows. In section [2] we briefly introduce the 
Higgs sector of the scale invariant NMSSM. In section[3]we present in detail the parameter values 
which we choose for the scan in the NMSSM parameter space. This is followed, in section [4] 
by the discussion of the SUSY particle effects in the loop mediated processes of the dominant 
NMSSM Higgs production through gluon fusion on the one hand and the decay into a photon 
pair on the other hand. Section [5] contains the numerical analysis with the presentation and 
discussion of X-k and mass distributions, of total widths, branching ratios and reduced rates 
with their correlations. A comparison with the present LHC Higgs search results is presented. 
Section [6] summarises and concludes the paper. 



2 The NMSSM 

We restrict ourselves to the NMSSM with a scale invariant superpotential. We do not take into 
account other possible extensions as the minimal non- minimal supersymmetric SM (MNSSM), 
new minimally-extended supersymmetric SM or nearly-minimal supersymmetric SM (nMSSM), 
neither extensions with additional U(l)' gauge symmetries 44], nor the case of explicit CP 



violation 138,45,46 



Including only the third generation fermions, the NMSSM superpotential in terms of (hatted) 
superfields is given by 

W = XSH u H d + ^S 3 + h t Q 3 H u t c R - h b Q 3 H d b c R - h T L 3 H d ? R . (2.3) 

The first term replaces the /i-term iiH u Ii d of the MSSM superpotential, while the second one, 
cubic in the singlet superfield, is introduced to break the Peccei-Quinn symmetry [47| in order to 
avoid the appearance of a massless axion. The last three terms represent the Yukawa interactions. 
The scalar mass parameters for the Higgs and sfermion scalar fields which contribute to the soft 
SUSY breaking Lagrangian read in terms of the fields corresponding to the complex scalar 
components of the superfields, 

-£mass = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + m||5| 2 

+ ™ 2 q 3 \QI\ +™$ R \tk\ + m ~ bR \~ b R\ +m%\Ll\ + r4 R \f 2 R \ . (2.4) 

And the trilinear soft SUSY breaking interactions between the sfermion and Higgs fields are, 

1 , ~ ~ 

- Aril = XA x H u H d S + — kA k S + h t A t Q 3 H u t c R - h b A b Q 3 H d b c R - h T A T L 3 H d f R + h.c. . (2.5) 

We work in the unconstrained NMSSM with non-universal soft terms at the GUT scale. The 
three SUSY breaking masses squared for H u , H d and S which appear in £ maS s can be expressed 
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through their VEVs by exploiting the three minimisation conditions of the scalar potential. 
While the MSSM Higgs sector at tree-level can be described by only two free parameters (in 
general chosen to be the mass of the pseudoscalar Higgs boson and tan /3) , the Higgs sector of 
the NMSSM is parameterised by the six parameters 

A , k , A x , A K , tan/3 = (H u )/(H d ) and fi e g = X(S) . (2.6) 

The brackets denote the VEV of the respective field inside. The sign conventions are chosen 
such that A and tan/3 are positive, whereas k, A\, A k and /x e g can have both signs. 

The Higgs sector consists of 3 CP-even Higgs bosons Hi (i = 1,2,3), two CP-odd states 
Aj (j = 1,2) and two charged Higgs scalars . The neutral Higgs bosons are ordered by 
ascending mass with H\ {A\) being the lightest CP-even (odd) Higgs boson. As higher order 
corrections to the Higgs sector are important and have to be considered in order to calculate 
the Higgs sector as accurately as possible, also the parameters from the non-Higgs sector, which 
enter through the loop corrections, have to be specified. These are the soft SUSY breaking mass 



terms in Eq. (2.4) for the scalars as well as the trilinear couplings in Eq. (2.5) and the gaugino 

(2.7) 



soft SUSY breaking mass parameters given by 

1 



gaugmos — ^ 



M 1 BB + M 2 ^W a W a + M 3 Y J G a G a + h.c. 



a=l a=l 



3 The Scan 

In the following we will perform a scan in the NMSSM parameter space in order to investigate 
the Higgs sector in view of the recent LHC Higgs search results together with the resulting 
possible theoretical and phenomenological implications. When performing our scan we seek to 
generate a Higgs spectrum where one of the scalar Higgs bosons corresponds to a state with mass 
value around 126 GeV leading to event rates in its production which are compatible with the 



LHC results. We furthermore keep the fine-tuning [18 19 26 27 as low as possible by demanding 
light top squark masses and/or small mixing in the stop sector. 

For the calculation of the SUSY particle and NMSSM Higgs boson spectrum and branching 



ratios we use the program package NMSSMTools [48 49 . The higher order corrections to the 
NMSSM Higgs boson masses are important [50| and have been included in NMSSMTools up 
to 0{atOL s + ctb a s) for vanishing external momentum. Within the package the Fortran code 



NMHDECAY 48 , an NMSSM extension of the Fortran code HDECAY 51,52 , provides the Higgs 



decay widths and branching ratios, while the SUSY particle branching ratios are obtained from 



the Fortran code NMSDECAY 1 53] based on the generalisation of the Fortran code SDECAY |52[|54 
to the NMSSM particle spectrum. The output of the NMSSM particle spectrum, mixing angles, 
decay widths and branching ratios is provided in the SUSY Les Houches Accord (SLHA) format 
[55] . Being interfaced with micrOMEGAs [56] , also the relic abundance of the lightest neutralino 
Xi as the NMSSM Dark Matter candidate can be evaluated with NMSSMTools. Furthermore, 
the package checks for the constraints from low-energy observables as well as from Tevatron and 
LEP. For details, we refer the reader to the program webpage |49|p| 

In order to restrict the parameter range for our scan we are guidelined by the following 
objectives which follow from theoretical and experimental considerations: 



1 Concerning the value of g — 2, it is non-trivial to find parameter combinations which can explain the 2cr 
deviation from the SM value. In our analysis we do not further consider this constraint. 



4 



• To keep the Higgs mass corrections (governed by the corrections from the (s)top sector) 
and hence the amount of fine-tuning as low as possible, the tree-level mass of the lightest 
Higgs boson is maximized by fixing tan /3 to small values chosen as 

tan/3 = 2, 4 . (3.8) 



Also the effective ^ e ff parameter is kept as low as possible in order to avoid fine-tuning. It 
is varied in the range 

100 GeV < //eg < 200 GeV . (3.9) 

Although we did not further consider the constraint coming from the anomalous magnetic 
moment of the muon, we decided to take positive values of // e ff as, similarly to the MSSM 



\i parameter, positive values are favoured when this constraint is included, see e.g. 11 



We shall be interested exclusively in large values of A in order to increase the tree-level 
mass of the CP-even Higgs boson associated with the 126 GeV Higgs boson resonance. 
At the same time we pay attention that it remains small enough to ensure the validity 
of perturbation theory up to large scales, chosen to be the GUT scale here. This also 
constrains possible values of n. Based on the results from the two-loop renormalisation 



group running down to 1 TeV with and without the possibility of exotic extra matter 19 
we hence perform our scan in the ranges 

0.55 < A < 0.8 and 10~ 4 < k < 0.4 . (3.10) 

The soft SUSY breaking trilinear couplings A\ and A K are varied in the ranges 

- 500 GeV < A K < GeV and 200 GeV < A x < 800 GeV . (3.11) 

For fine-tuning reasons we keep the soft SUSY breaking masses of the stop sector rather 
low and vary them simultaneously as 

500 GeV < Mq 3 = M- tR < 800 GeV . (3.12) 

For Ay (U = u,c, tj^we choose two representative values corresponding to low and large 
mixing, 

Au = GeV and 1 TeV . (3.13) 
Our lightest stop mass is hence about 400 GeV and in accordance with the LHC constraints 



57 



• In order to comply with the present LHC search bounds 60 1 , we conservatively set the soft 
SUSY breaking masses of the squark sector of the first two generations equal to 2.5 TeV 
and, for simplicity, also those of the slepton sector apart from the soft SUSY breaking stau 
masses. The latter are chosen equal to 300 GeV. This way we still allow for rather light stau 

2 In NMSSMTools there is no distinction between A u , A c , At. 

3 In scenarios with a very small mass difference between the lightest stop ii and the lightest neutralino Xi 
assumed to be the lightest SUSY particle, stop masses down to about 100-130 GeV are still allowed for m^o > 
90 GeV [58l|59l. 
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masses but are conservative enough to fulfill the latest LHC results 61 . It should be noted, 
however, that our results almost do not change by choosing different values in the stau 
sectoij^] as the influence of the slepton sector on the Higgs mass corrections is negligible. 
And contrary to the MSSM, light stau masses here do not lead to an enhancement of 
the partial width into photons [fjj, as we have chosen small values of tan/3 and fi e g. We 
furthermore set the trilinear couplings of the down and lepton sector equal to 1 TeV and 
the right-handed soft SUSY breaking sbottom mass equal to 2.5 TeV. This results in 
light sbottom masses of about 500 GeV < < 800 GeV. Hence we have (D = d, s, b, 
E = e,n,r) 

Mu R = M- CR = M 3r = Mq i 2 = M iR = M~ m = M~ Li 2 = 2.5 TeV, 

M fR = M Lz = 300 GeV , A D = A E = 1 TeV . (3.14) 

• The gluino soft SUSY breaking mass parameter has been set to 

M 3 = 1 TeV . (3.15) 

The remaining two soft SUSY breaking gaugino parameters have been chosen M\ = 
150 GeV and M 2 = 300 GeV. 

It should be noted that in NMSSMTools the NMSSM-specific input parameters A, k, A\ and A K 
according to the SLHA format are understood as running DR parameters taken at the SUSY 
scale M = 1 TeV, while tan/3 is taken at the mass of the Z boson, Mz- 

We remark, that at the cost of a more time consuming scan we could have enlarged our 
parameter ranges of A K , A\ and k. As will be evident from our numerical analysis later, the 
limitation of the scan to this restricted parameter area nevertheless leads to a substantial amount 
of parameter points which are compatible with the applied constraints due to experimental 
results and fine-tuning arguments. Note also, that choosing large positive values for A K for 
negative k leads to non self-consistent solutions. Concerning A\, it is related to the charged 
Higgs boson mass, which is below the experimental limit if A\ is chosen too small. A posteriori 
it also turned out that the chosen upper bound of A\ was largely sufficient to capture the 
maximum of allowed parameter points which can be achieved for the chosen A K range. 

The parameter scan is further restricted by demanding the NMSSM Higgs spectrum to fulfill 
the following conditions: 

• We demand one of the scalar Higgs bosons, which we will denote from here on by h, to 
have its mass in the range 

scalar Higgs boson h: 124 GeV < m h < 127 GeV , (3.16) 

where we have conservatively assumed a 3a error on the mass value of the scalar particle 
discovered at the LHC fiM. 



In order to explore the possibility of an enhanced branching ratio into photons, we fur- 
thermore demand that the 77 rate around the invariant mass value 126 GeV fulfills: 

rate for the 77 final state normalised to the SM value > 0.8 . (3-17) 



4 Also the SUSY breaking masses of the squarks of the first two generations barely influence the outcome of 
the scans. 
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We do not put any restrictions on the rates in the massive gauge boson and fermion final 
states. 

For the other Higgs bosons, i.e. the pseudoscalar Higgs bosons and the scalar Higgs bosons 
outside the mass range around 126 GeV, we check if they have not been excluded by the 
LEP, Tevatron and LHC searches. Otherwise the whole parameter point is rejected. We 
have taken into account the newest exclusion limits in the various final states reported by 
the experiments [62|470| , which we have implemented in NMSSMTools. 



Our choice of parameters results in rather low top squark, charged Higgs boson and chargino 
masses, still compatible, however, with present LHC SUSY search results. For tan (3 = 2 we 
have 

mfi = 400 - 820 GeV , mg, = 530 - 890 GeV , (3.18) 
M H ± = 200 - 500 GeV , M~± = 105 - 165 GeV , M~± = 345 - 360 GeV , (3.19) 

and similar values for tan f3 = 4. The stop mass values are small enough so that the fine-tuning 
is expected to be rather low. 

We finally remark that we did not restrict our parameter points taking into account the relic 
density. We checked, however, that there is a substantial amount of parameter points which lead 
to relic densities due to a neutralino DM candidate, which are smaller than the WMAP value. To 
achieve the correct amount of relic density another candidate than the neutralino would have to 
be thought of. Furthermore, we convinced ourselves that e.g. by slightly changing the values of 
the gaugino mass parameters Mi, M2 the correct amount of relic density could be achieved, while 
the Higgs mass spectrum remains practically unchanged, so that we did not further consider 



this constraint. For discussions taking into account DM constraints, see e.g. 21-24,37 . 



4 NMSSM Higgs boson production and decay 

In order to decide whether the 126 GeV NMSSM Higgs boson reproduces the rates as measured 
by the experiments, its production cross sections and branching ratios have to be investigated. 
In the following the dominant production process through gluon fusion and its modification with 
respect to the SM will be discussed in detail. We furthermore investigate the NMSSM Higgs 
branching ratio into photons, as the LHC experiments see a slight excess here with respect to 
the SM. With the presently available data, this has to be taken with due caution, however, as 
it could still turn out to be a statistical fluctuation. If it persists, however, it is a hint towards 
New Physics and shall be taken into account in our analysis. We start with some preliminary 
remarks and set up our notation. 

At the LHC, for small values of tan /3, the production processes for a single neutral CP-even 
NMSSM Higgs boson Hi (i = 1, 2, 3) or a CP-odd Higgs state Aj (j = 1, 2) are given by 

Gluon fusion: gg — > Hi and gg — > Aj 

Gauge boson fusion: qq — >■ qq + W*W* /Z*Z* — > qqHi 

Higgs-strahlung: qq Z* /W* ->■ Hi + Z/W 

Associated production with tt: gg/qq — > ttHi and gg/qq — )• tiAj 
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where gluon fusion is the most important process followed by gauge boson fusion. Higgs- 
strahlung and associated production with a top quark paiij^] only play a minor role and are 
more important for the determination of Higgs boson couplings. 

The NMSSM production processes and decay channels deviate from the corresponding SM 
Higgs H SM processes due to modified Higgs couplings and additional SUSY particles running in 
the loop mediated processes. The couplings of the CP-even Higgs states Hi (and also those of 
the CP-odd states) depend on their decompositions into the weak eigenstates Hd, H u and S, 

Hi = Si t d Hd + S\ :U H u + Si iS S , 

H2 = S24 Ha + S2, u H u + S2, s S , (4-21) 
Hs = S34 Hd + S^u H u + Ss jS S . 

The coefficients Si )U , Si : d hence quantify the amount of up- and down-likeness, respectively, while 
Si S is a measure for the singlet-component of a Higgs mass eigenstate. Mixings between the 
££/(2)-doublet and singlet sectors are proportional to A, and can be sizeable for A > 0.3, leading 



to significant effects on the NMSSM Higgs couplings and hence phenomenology 15, 16 18, 19 . 

The inclusive production cross section o"j nc ; for a CP-even Higgs boson is composed of gluon 
fusion, vector boson fusion, Higgs-strahlung and associated production with ti, 

a ind (H) = a(gg -> H) + a(Hqq) + a{WH) + a(ZH) + a(tiH) « a(gg -> H) , (4.22) 

with H = Hi,H SM , respectively. It is dominated by the gluon fusion cross section. For later 
convenience in the discussion of our results we normalise the relevant quantities of the NMSSM 
Higgs bosons to the corresponding SM counterparts. Thus we define the ratio Ra incl of the 
NMSSM inclusive cross section compared to the SM one, 

Ra mcl (Hi) = ^J^SM) ~ RvggiHi) ' ( 423 ) 

where we have used R agg {Hi) defined as the ratio of the NMSSM gluon fusion production cross 
section to the SM one, 

If not stated otherwise, in these and the following ratios the mass of the NMSSM Higgs boson 
Hi and the one of the SM Higgs H SM are taken to be the same and they are subject to the 
constraint M hSM = M Hi =m h = 124 - 127 GeV. 

The ratio Rr tot for the total width compared to the SM Higgs total width is given by 

Rr tot m = ^§iry (4-25) 

While in the SM the largest decay width of a Higgs boson of about 126 GeV is the one into bb, 
the most important search channels are given by the 77, the massive gauge boson and the rr 
final states. We define the ratios of the NMSSM Higgs decay partial widths relative to the SM 
as {X = 7, W, Z, b, t) 



3 For small tan /3 values associated production with a bottom quark pair is negligibly small. 
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The ratios of branching ratios are given by 

r>BR ( tt \ — BR(Hj ->■ XX) _ R Txx {Hj) 

K xx^i) = BR{H su xxj ~ R Vtotm ■ (4 - 27j 

The experimentally observed rate in a given channel X is given by the reduced cross section Rxx 
which is obtained from multiplying the Higgs production ratio relative to the SM, R a . ncl (Hi), 
with the Higgs branching ratio for the channel of interest relative to the SM. For example, for 
the two photon final state we have 

R^m = Rcr incl (Hi)R^ (Hi). (4.28) 

The corresponding reduced cross sections in the other decay channels VV (V = W,Z), bb, tt 
may be similarly expressed, namely: 

Rvv(Hi) = R atncl (Hi) Ryy(Hi), R^(Hi) = R a (VH)(Hi) R^ b R (Hi), 

R TT m = R^mR^Hi). (4.29) 

In the bb final state we restrict ourselves to associated production of the Higgs boson with a 
W or Z boson, as we will compare our results later with values given for this channel by the 
experiments. 

It is important to note that there can be NMSSM spectra where two neutral Higgs bosons lie 
close in mass. Due to the limited experimental resolution these cannot be separated from each 
other and both contribute to the signal. The program NMSSMTools takes this into account by 
super-imposing the signal from the nearby Higgs boson with a Gaussian weighting. The width 
of the Gaussian smearing is adapted to the respective experimental resolution in the different 
final states, where clearly the 77 and ZZ final states have the best resolution, while the mass 
resolution in the tt and bb final states is less good, and in the WW final state the mass cannot 
be reconstructed. Hence, the ratios for the rates, Rxx, depending on the scenario and related 
NMSSM spectrum under consideration, can be superpositions of rates of different Higgs bosons. 
In favour of an unambiguous notation and to make contact with the signal strengths \i = a j ctsm 



reported by the LHC experiments, we denote by \ixx the reduced cross sections (4.28), (4.29), 
which are built up by the superposition of the rates from the 126 GeV h boson and another 
Higgs boson <3? = Hi,Aj, which is close by in mass, 

lixxih) = R a {h) R X x(h) + Yl R*@)Rxx(*)F(M h ,M*,d XX ) ■ (4-30) 

* ^ h 
\M*-M h \ < S 

Here a = a(VH) in case X = b and a = Oi nc i otherwise. By 5 we denote the mass resolution 
in the respective XX final state and by F(Mf l ,Mj„dxx) the Gaussian weighting function as 
implemented in NMSSMTools. The experimental resolution of the different channels is taken into 
account by the parameter dxx, which influences the width of the weighting function. We impose 



the restriction Eq. (3.17) on the thus calculated 77 rate, which in fact is the one observed in 



experiment. Hence, in summary the conditions we impose on our parameter points are: 
Conditions on the parameter scan: 

At least one CP-even Higgs boson h with: 124 GeV < M h < 127 GeV 
The reduced cross section for 77 must fulfill: /U 77 (/i) > 0.8 with (4-31) 

124 GeV <M h = M hSM < 127 GeV 
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4.1 Higgs boson production through gluon fusion 

The cross section for NMSSM Higgs production via gluon fusion is mediated by quark Q and 
squark Q triangle loops, cf. Fig. [l| The latter become particularly important for squark masses 



below about 400 GeV |4[|5[|9[|7l|. At leading order (LO) in the narrow- width approximation the 
hadronic cross section for scalar Higgs bosons Hi (i = 1, 2, 3) can be cast into the form [7l|-[73] 



dr Hi 



cr, 



Hi 



G F a 2 s (n R ) 
288\Z2vr 



(4.32) 
(4.33) 



with the gluon luminosity C gg , the Fermi constant Gf, Tu i = M^./s, where s denotes the 
squared hadronic cm. energy and Tx = 4:M X /Mjj. (X = Q,Q). The strong coupling constant 
a s is taken at the scale chosen equal to the mass of Hi. The form factors A 1 ~ are given by 



-Ml 



;i - r)/(r)] 



:r[l-rf(r)] 



(4.34) 
(4.35) 



and the function f(r) reads 



f(r) 



1 



arcsm 



log 



r 

1 + VT 



T 



ITT 



T > 1 



T < 1 . 



For large values of the loop particle masses the form factors become constant, 

r2 ^ h™2 



1 
1 

4 



for Mjj. < 4mg 
for Mtf. < 4m^ . 



(4.36) 



(4.37) 
(4.38) 



For small values of tan f3 the most important contributions come from the top and stop loops. In 
order to study the effect of the stop loops and their interplay with the top quark loop, the Higgs 
couplings to the top and stop quarks, gq^g^ 1 , have to be investigated. Due to the diagonal 
gluon coupling to stops, in the loop only the Higgs couplings to two equal stops can appear. 




9 anr> 



Hi 



Q' > Hf 

I / 




Hi 



Figure 1: Generic diagrams contributing gluon fusion production of Hi 
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Hence for g-* we have to consider the couplings 



(Si * cos P - S{ u sin /3) ( ~~ cos 2 Q~ t - \ sin 2 6 W cos 2%) + 

m~ 2 3 sinn 

ti ti 



At; cos /3 

/*eff"i,d "t" -A-tDi^u ,~>, 



1,3 



1 m* 
+-sin2fl f ~ 

2 l m sm/3 

M 2 1 2 r7i 2 9- 

(5 M cos /3 - 5 JjU sin /3)^(- sin 2 t - + - sin 2 % cos 2%) + * .'' M 

12 



(4.39) 



m? 2 

*2 



1 n m t 

- sm2%— o — ; — r 

2 r mJ sm/3 

12 



At; cos P 



(4.40) 



The Higgs mixing matrix elements Si jX (x = d,u,s) have been defined in Eq. (4.21). We are 



interested in NMSSM parameter scenarios with one of the CP-even Higgs bosons having a mass 
around 126 GeV and production rates which are not too far away from the corresponding SM 
rates in the various final states in order to comply with the LHC Higgs search results. For the 
sake of simplicity in the investigation of the couplings we therefore define the SM limit of the 
NMSSM, which is given by first performing the MSSM limit, which is recovered by A, k — > 
with k/A constant and keeping the parameters fi e s,A\ and A R fixed. Within the MSSM limit 
then the decoupling limit is performed. In the thus defined SM limit the mixing matrix elements 
of h become 



cos/3 



sin p 







(4.41) 



And we get for the couplings to the stops 



Ml A 



m~ "2 
ti 



cos 2/3 ^(^ cos 2 6 i - ^ sin 2 6 W cos 2%) + ^ 



ra- 
ti 



m t 



+-sin20r- „ 
2 l mf sin/3 



-/x e ff cos P + A t sin /3] 



g ht 2 t 2 



M 2 1 2 to 2 

cos 2p—f(- sin 2 Of + - sin 2 8 W cos 29 1) H 1 

m~ 2 3 m 



sin 2Q- t „ m f 

m| sm p 
t 2 



-/x e ff cos P + A t sin /3] 



(4.42) 



(4.43) 



In the scenarios of the parameter scan which are left over after applying our criteria Eq. (4.31 ), 
for tan/3 = 2 the relations (4.41) are approximately fulfilled apart from the singlet component, 
which can take values of up to ~ 0.1. The approximation gets worse in scenarios with strong 
singlet-doublet mixing, where the singlet component can reach values of up to ~ 0.6. In this case 
we can have suppressed couplings of the 126 GeV Higgs boson to bottom quarks. For tan/3 = 4 
the behaviour is similar for large mixing, while for small mixing the deviations from this SM 
limit are more important. Note, that suppressed couplings to top quarks are largely ruled out 
due to our demand of the 77 reduced cross section exceeding 80% of the SM value. This can 
only be achieved if the dominant production cross section through gluon fusion is large enough, 
which is not the case for Higgs couplings to the top quarks being too suppressed compared to 
the SM value. Nevertheless, also the top quark couplings can be suppressed compared to the 
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SM in the cases where the branching ratio into 77 is enhanced or where the Higgs rates are 
built up by the superposition of rates stemming from more than one Higgs boson, so that our 
restriction on the 77 rate can be fulfilled. 

The mixing angle 6f which diagonalises the stop mass matrix is given by 



sin 2% 



2m t {A t - /j, eS / tan/3) 



m 



ti 



in - 



(4.44) 



*2 



where m j 1(2) denotes the lighter (heavier) stop quark mass. For the sake of the discussion we 



assume fi e g to be zeroj^] With this approximation and neglecting small D-term contributions we 

2mtAt- For large values of At, the m 2 



have for M 



Qs 



Mr the mass difference m~ 



2 
ti 



term and the last term in (4.42) therefore have opposite sign (in (4.43) same sign). Neglecting 



the small D-texm contribution given by the first term in the coupling, the Higgs coupling to the 
lighter stops t\ becomes negative in this case. Assuming t\ to be relatively light, the contribution 
from the t\ loop will be more important and we will not consider the ti loop contribution in 
this case. The Hi coupling to the top quarks on the other hand is given by 



Hi 
9Q 



9Hdt 



sin/3 



which becomes in the SM limit for the 126 GeV Higgs boson 

9htt = 1 • 



(4.45) 



(4.46) 



Hence for large values of At the t\ and the top loop contributions interfere destructively so 
that the gluon fusion cross section decreases. For small values of At there is no mixing in 
the stop sector leading to a positive Higgs coupling to t\ and constructive interference, thus 
enhancing the gluon fusion production cross section. For non-zero intermediate At values, the 



last two contributions of (4.42) cancel each other (the exact value of At depends on the specific 
parameter choices), and the stop contribution is smallQ The Higgs coupling values to the 
sbottoms are hardly influenced by a change in At which enters in the mixing matrix elements 
Si )X only through higher order corrections to the Higgs boson masses. 



4.2 Higgs decay width into two photons 

The decays of the scalar NMSSM Higgs bosons into photons are mediated by W boson and heavy 
fermion loops as in the Standard Model and, in addition, by charged Higgs boson, sfermion and 
chargino loops; the relevant diagrams are shown in Fig. [2] The partial decay widths, adapted 
from the MSSM result (7l}[73}[74) , are given by 



T(Hi -> 77) 



G F a 2 M%. 
128^/2vr 3 



J2N cf e 2 f gfAf(T f )+g§A%(T W ) 



+ 9% ± A h h %(th±) + £ <$4?i(r i± ) + Y d N cf e)gfAf{r } ) 



, (4.47) 



6 Our values of /x e g- are small enough not to change the conclusions for non-zero values. 

7 The influence of the stop loop contributions on gluon fusion and decay into photon final states has been 
discussed in the context of the MSSM in detail in jHJ. 
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Figure 2: Generic diagrams contributing to the decay iifi — >• 77. 



with the colour factor N c f = 1(3) for leptons (quarks) and et denoting the electric charge of the 
loop particle. The form factors are given by 



< l (r) 



2r[l + (l-r)/(r)] 

-r[l-r/(r)] 

-[2 + 3r + 3r(2-r)/(r) 



(4.48) 
(4.49) 
(4.50) 



where r = 4M|/M|. with being the mass of the particle X in the loop. For large loop 
particle masses Mj the form factors approach constant values, 



A h^ 



-7 



for Ml. <C 4M 2 _ ± 
for Ml. < 4M 2 

1 H ± ,f 

for M|. < 4M^ . 



(4.51) 



The Higgs couplings to fermions, W bosons, charged Higgs bosons and charginos, appearing in 
the decay width into two photons, are given by 



Hi 

Hi 
9w 



Hi 
9 H ± 



Si^ u / sin/3 for / = up-type fermion 
•S^d/ cos/3 for /= down- type fermion 

Si a cos /3 + Si u sin/3 



Hi H{ 

9 * ±=g xtxZ 



2M W 
~M^ 



_qkkSi t d + 



&^(cos 3 /3 5 M + sin 3 /3^, u ) 
)s z 8w V / 



cos 

2 cos 2 9w 



+ - cos (3 sin /3( (3 + tan 2 6*iy) — 4A 2 /g 2 J ( sin/3 S 1 ^ + cos/3 iSj )U ) 



(2A/i cff + sin 2(3(A X \ + 2^ eff )) S, 



(4.52) 

(4.53) 
(4.54) 



(4.55) 



The matrix elements qj~i, Ski,fki (k, I = 1,2) in terms of the matrix elements of the matrices U, V 
diagonalising the chargino mass matrix [75] read 



Qkl 



1 

71 



U a V k 



1 j 



1 

71 



Aw 



(4.56) 



H * 
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In the SM limit, denned in the previous subsection, the couplings become 



Hi 

Hi 
9 H ± 



Hi 

9-i-± 

xtxt 



1 
1 



C ° S(2 9 ^ ) ( cos 4 P + sin 4 ^ + cos 2 /3 sin 2 /3 f (3 + tan 2 ^) - 4- 
2 cos^ Bw V / V g 



2MW ( R^ ■ R 

{q kk cos (3 + s kk sin /3 



(4.57) 
(4.58) 

(4.59) 
(4.60) 



The sfermion loop contributions are only important for the third generation with light masses. 
While it has been shown in the context of the MSSM that the f loop contributions can enhance 
the decay width into photons [6j , this is not the case here as we assume small values of tan /3 and 
fi e ff. Due to the small tan/3 values also sbottom loops play a minor role. The most important 
sfermion contribution comes from the stop loops in particular a light t\. The Higgs coupling to 



the stops has been given in section 4.1 Eqs. (4.39), (4.40) and for the SM limit in Eqs. (4.42) 



(4.43). Since in the Higgs decay into photons the quark and W loop contributions interfere 



destructively, in the decay the effect of the stop loops is opposite to the one in the production. 
For At = GeV the stop loop contribution suppresses the decay into photons, for A t = 1 TeV 
it leads to an enhancement. 

For light enough chargino and charged Higgs boson masses their loop contribution also plays 
a role. We find that they can lead to an enhancement for the partial decay width into photons 
(see also fflpjpl). 



5 Numerical Analysis 

In this section we show our numerical results. When performing the scans we find scenarios 
in which both the lightest scalar Higgs H\ and the heavier one H2 can have masses around 
126 GeV. We will call the respective Higgs boson in this case h. Not all its couplings are 
necessarily equal or near the corresponding SM Higgs coupling values. More importantly the 
reduced cross sections in the various final states, which can be superpositions of signals from 



Higgs bosons with masses close by, as defined in Eq. (4.30), have to reproduce the experimental 
results. To avoid a flood of plots, in the following we will only show the ones for tan/3 = 2 and 
comment on the plots for tan /3 = 4. 



5.1 Parameter values and mass distributions 

In Fig. [3] we show the distributions of the allowed parameter points in the k—X plane leading 
to Hi representing the CP-even Higgs boson with mass in the range 124-127 GeV (left) and to 
H2 (right) being h, respectively, for the two values of At = GeV and 1 TeV. The colour code 
denotes the number of points. As can be inferred from the figures, we have much more points 
allowing for h = H2 than for h = H\. This can be explained as follows. As we demand the stop 
mass parameters to be rather low, which leads to smaller higher order corrections for a fixed 
mixing, the parameter A at low energies must be large enough to get to the right mass. The 
demand of perturbativity up to the GUT scale then implies stringent constraints on the coupling 
ft(l TeV). In addition /u e g is required to be smaller than 200 GeV to avoid fine-tuning. Because 
the masses of the extra NMSSM scalar and pseudoscalar states, which are predominantly SM 
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Figure 3: The distribution of allowed parameter points in the k-A plane for h = Hi (left) and h = H2 (right) for 
A t = GeV (upper) and At = 1 TeV (lower). The red (blue) contour lines show the two-loop upper bounds for 
A and k at 1 TeV in the NMSSM without (with) extra matter above 1 TeV. The colour code denotes the number 
of points. 



singlets, are set by Kfj, e ff/\ these states tend to be lighter than the 126 GeV Higgs boson. As 
both the Hi and H2 mass values increase with rising k, for Hi we need for the same reasons 
large k values of k ~ 0.4, while too large values of k lead to too large H2 masses so that here 
values k ~ 0.07-0.09 are preferred. With increasing values of At the stop mass corrections to 
the tree-level masses become more important so that a 126 GeV Higgs mass can be attained 
more easily and therefore more parameter points pass the constraints. For the same reason the 
maximum of points is given for smaller values of A now, decreasing from A ~ 0.73 (0.72) at small 
stop mixing to A ~ 0.68 (0.66) for h = Hi (H2) at large mixing. 

In the plots we also show the upper bounds on A and k imposed by perturbativity derived 
from the two-loop renormalisation group running from the GUT scale down to 1 TeV. These 
limits can be somewhat relaxed when allowing for extra exotic matter with mass around 1 TeV. 
They show that an Hi Higgs boson with mass around 126 GeV can only be achieved for large 
mixing with At = 1 TeV. For lower values of At even with the inclusion of extra matter, this 
is not possible. The heavier Higgs boson H2 on the other hand can have a 126 GeV mass 
value with and without exotic matter. We finally note that in case At = GeV, for h = Hi 
the trilinear couplings A K , A\ cluster around (A K ,A\) = (0 GeV, 310 GeV) and for h = H2 
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Figure 4: Mass spectrum of Ai and H2 for h = Hi (left) and of Ai and Hi for /i = H2 (right) with At = GeV 
(upper) and 1 TeV (lower). The colour code denotes the number of points. 



around (A K ,A X ) = (-140 GeV, 310 GeV). In case A t = 1 TeV, we have for h = Hi the 
maximum of points around (A K ,A\) = (0 GeV, 340 GeV) and for h = H2 around (A K ,A\) = 
(-140 GeV, 340 GeV). 

Figure [4] shows the mass distributions of the lighter neutral Higgs bosons for Hi and H2 
being h, respectively. For h = Hi there exist parameter regions where Hi and/or Ai are very 
close in mass. Depending on the respective experimental resolution in the investigated final state 
their signal can superimpose the h rate. This superposition has been taken into account in the 
reduced cross sections discussed later. The maximum of parameter points clusters around mass 
values Mh 2 ,Ai ~ (175, 170) GeV. Also for Hi with mass ~ 126 GeV the Hi and/or A\ state can 
be close in mass and contribute to the signal. Their masses can be also much smaller, however, 



so that Hi decays into these final states become possible, leading to distinct signatures 76 . The 



maximum parameter points are found for Mh x Ai ~ (85, 110) GeV. The masses of the heavier 
Higgs bosons H3 and Ai lie between about 300 and 500 GeV. 

We remind the reader that in all plots we have already taken into account the latest exclusion 
limits from LEP, Tevatron and LHC which apply to the non-/i Higgs bosons. In particular for 
scenarios with h = Hi this leads to a substantial reduction of allowed parameter points. 

As for tan/3 = 4, it turns out that for small mixing no parameter combination fulfills the 
conditions (4.31 ) for the lightest NMSSM Higgs boson Hi. Only for large mixing a few hundred 
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Figure 5: Branching ratio into 77 relative to the SM against the normalised total width for h = Hi (left) and 
h — H2 (right) and At = 1 TeV. The orange points (for h = H2) indicate that non-SM decays are allowed. 



parameter points survive which cluster around (k, A) = (0.4,0.8). For h = H2 both in the 
low and in the large mixing case the conditions are fulfilled with the maximum of points in a 
somewhat extended region around (re, A) = (0.4,0.8). The reason is that for larger values of 
tan/3 the tree- level upper mass bound is lower than for tan/3 = 2, so that more substantial 
higher-order mass corrections are needed which in case H\ is to have a mass around 126 GeV 
can be achieved only for large mixing. If we now apply the perturbativity bounds on k and A 
it turns out, however, that none of the h = Hi scenarios survives as even with extra matter at 
1 TeV the maximum allowed value is A = 0.66 for k = 0.4. For the h = H2 scenarios a few 
scenarios survive if extra matter is included. Otherwise the perturbativity bounds imposing a 
maximum value A = 0.64 for k = 0.4 are not respected. 

We summarise, that for tan/3 = 2 there are scenarios with h = H\ which respect perturba- 
tivity in case of large mixing and inclusion of extra matter at 1 TeV. For h = H2 this is the case 
for both small and large mixing and NMSSM with and without extra matter. For tan /3 = 4 
only h = H2 scenarios survive and are compatible with perturbativity for low and large mixing 
if extra matter is included. 

5.2 The 77 final state 

We first discuss the behaviour of the photonic branching ratio, which is shown in Fig. [5] with 
respect to the SM compared to the normalised total width for h = H\ and H2, respectively, and 
At = 1 TeV. The plots for At = GeV look very similar and therefore are not shown here. The 
branching ratio into 77 can be largely enhanced up to ~ 5.6 times the SM value. This is due to a 
substantially suppressed total width because of strong singlet-doublet mixing of the Higgs boson 
with mass ~ 126 GeV. Its coupling to bottom quarks is therefore strongly reduced, leading to 
a small total width (dominated by the decay into 6-quarks) and hence an enhanced branching 
ratio. The increase in the branching ratio, however, can also be due to an enhanced decay width 
into photons caused by squark, charged Higgs and/or chargino loop contributions, as has been 
discussed above. Therefore also for enhanced total widths the branching ratio can be larger 
than in the SM case. If, however, besides the couplings to bottom quarks also the other Higgs 
couplings are substantially suppressed due to strong singlet-doublet mixing, the loop-induced 
coupling to photons becomes very small, leading to small branching ratios also in the case of 
small total widths. For h = H2 we can observe that the total width can be increased by up 
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Figure 6: Branching ratio into 77 relative to the SM against the inclusive production cross section relative to the 
SM for h = Hi (left) and h = H 2 (right) for A t = GeV (upper) and At = 1 TeV (lower). For green/red points 
the perturbation theory is valid up to the GUT scale. Cyan/pink points require extra matter above 1 TeV and 
yellow points violate the two-loop upper bounds on A, k both with and without extra matter. Points above the 
black line lead to an enhanced R~,~,. 



to ~ 1.5 compared to the SM. This happens where decays of H2 into other lighter Higgs bosons 
Hi or A\ and/or neutralino final states are kinematically allowed 76 . The relevant decays are 
H2 — > X1X11 H2 — > H\H\ and H2 — > A\A\, with the latter being rarely realised. 

The corresponding plots to Fig. [5] for tan (3 = 4 show a similar behaviour with altogether less 
parameter points, however, and a maximum photonic branching ratio enhancement of i2 77 ~ 5 
for both h = H2 and Hi (with only the large mixing case surviving here). And for the total width 
the maximum value is i?p. . ~ 1.35 due to H2 decays into light Higgs bosons or neutralinos. 

With Fig. [6] we discuss the interplay of production and decay on the photon rate. We show the 
branching ratio into 77 relative to the SM plotted against the inclusive cross section normalised 
to the SM for either h = Hi or h = H2. As the inclusive production is dominated by gluon fusion, 
we can restrict our discussion to this production process. The figures show that for vanishing 
At gluon fusion can indeed be enhanced compared to the SM due to stop loop contributions, 



as has been discussed in Section 4.1 With rising mixing the stop loop contribution interferes 
destructively, and for A t = 1 TeV the gluon fusion process is suppressed compared to the SM. 
Also the branching ratio into photons shows the expected opposite behaviour. For large values 
of At, where constructively interfering stop loops enhance the partial width, we can observe 
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Figure 7: Reduced cross sections into 77 against Mh ± (left) and M^o (right) for h — H2, tan /3 
The colour code denotes the size of the branching ratio BR(H2 — > H\H\) (left) and BR(H2 



2 and A t = 1 TeV. 
* XiX?) (right). 



slightly larger branching ratios than for At = GeV. It should be kept in mind though that the 
behaviour of the branching ratio is an interplay of the partial width into photons and the total 
width. Once again for tan/3 = 4 the corresponding plots to Fig. [6] show a similar behaviour with 
altogether less parameters points. 

Above the black line the reduced cross section R~~ = R^R a . , > 1 8 As can be inferred 

I I 7 / u incl — 1 I 

from the plots, in the NMSSM both Hi and H2 are compatible with a 126 GeV Higgs boson and 
an enhanced rate into photon final states. For h = H2 there are substantially more (red) points, 
which are compatible with the constraints that come from the requirement of the validity of 
the perturbation theory up to the GUT scale, than for h = H\ (green points). In particular for 
vanishing At extra matter is required, behaviour which can be traced back to the need of the 
H\ tree-level mass being as large as possible, cf. the discussion in the previous subsection. We 
note that there are scenarios where both the branching ratio and the inclusive production are 



very small due to h being very singlet-like. These scenarios passed the constraint (4.31) as in 
this case the photon reduced cross section /x 77 , which can be a superposition of contributions 
from various Higgs bosons being close in mass, is dominated by the contribution from another 
light Higgs boson with a mass of ~ 126 GeV, which is not singletdike in this case. 

As already mentioned above the heavier CP-even Higgs boson H2 can decay into a pair of 
lighter Higgs bosons or neutralinos in certain parameter regions. This is shown in Fig. [7j where 
for tan (3 = 2 and At = 1 TeV the reduced cross section in the 77 final state in case of h = H2 , 
^ 77 (i?2), is plotted against the mass of the lightest scalar Higgs boson Hi and the mass of the 
lightest neutralino Xi > respectively. The colour code denotes the size of the respective branching 
ratio, which is zero above the kinematic thresholds. These rainbow plots show that in case of 
enhanced photonic rates such non-standard Higgs decays always remain below about 10-20%. 
The reduced cross section /i 77 is suppressed in case of sizeable branching ratios above ~ 0.25 
with a maximum of BR^ x (HiHi) w 0.36 and BR^ x (xiXi) ~ 0.43. They are small enough 
not to be excluded by the present experimental bounds. As can be read off Fig. [7] (left), the 
largest enhancements in the photon final state occur for almost degenerate Hi and H2 masses, 
which corresponds to neutralino masses around 73 GeV, see Fig. [7] (right). We explicitly verified 

8 Note that we discuss here the reduced cross section for h only. Later we will look at reduced cross sections 
fJ-xx in the final state X, built up by the 126 GeV Higgs boson and possibly nearby Higgs resonances. This is 
what actually is observed in the experiment. 
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Figure 8: Reduced cross section into 77 versus the reduced cross section into ZZ for At = GeV (upper) and 
A t = 1 TeV (lower ) with h — Hi (left) and h = H2 (right). Cyan/pink points indicate the signals where at least 
two Higgs bosons with similar masses overlap and the combined reduced cross section deviates by more than 10% 
from the reduced cross section of the individual Higgs boson. Bars: Experimentally measured values with error 
bars (full/ ATLAS, dashed/CMS). 



that here the enhanced rate in the photon final state is due to the increased branching ratio into 
photons because of suppressed H2 couplings to b quarks in this case. Due to sum rules the Hi 
coupling to b quarks is then substantial. The combination of the effects of Higgs couplings to 
SM particles and experimental exclusion limits then implies the observed pattern in the plots. 
Concerning H%, it mainly decays into 6-quark pairs with a branching ratio of 0.8-0.9, followed 
by decays into r pairs and a branching ratio of roughly 0.1. The Higgs-to-Higgs or Higgs-to- 
neutralino decays hence lead to interesting final state signatures with e.g. 46, 2b2r, 4r or even 
multi-/i final states in the former case, from the secondary Higgs decays. In the latter case the 
final state lightest neutralino entails large missing energy. Such events could act as smoking gun 
signatures for extended Higgs sectors beyond the minimal SUSY version. 



5.3 Compatibility with the LHC Higgs search results 

In this subsection we investigate the compatibility of the results for the reduced cross sections 
fixx with the experimental best fit values of the signal strengths in the various final states. 
Figures [8 11 show the reduced cross section in the 77 final state compared to the one in ZZ, 
WW, bb and rr, respectively, for h = H x and H 2 with A t = GeV and A t = 1 TeV. The 
bars represent the newest results for the best fit values of the signal strengths fi = cr/asM 
in the different final states, reported by the ATLAS [lj |64}|66] and the CMS Collaboration 
(2 67-70 , together with their corresponding errors. The values and errors are listed in Table [T] 
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Figure 9: Reduced cross section into 77 versus the reduced cross section into WW for A t = GeV (upper) and 
A t — 1 TeV (lower) with ft = Hi (left) and ft = H2 (right). Cyan/pink points indicate the signals where at least 
two Higgs bosons with similar masses overlap and the combined reduced cross section deviates by more than 10% 
from the reduced cross section of the individual Higgs boson. Bars: Experimentally measured values with error 
bars (full/ ATLAS, dashed/CMS). 



in Appendix [Aj First of all the plots demonstrate that both Hi and H2 can have a mass around 
126 GeV and be compatible with the experiment, for small and for large mixing in the stop 
sector. Moreover an enhancement in the photon rate by up to a factor ~ 2.4 is possible. The 
allowed parameter regions are somewhat more extended for At = GeV, which is an interplay 
between the production cross section and decay into photon pairs leading to more important 
reduced rates for the small mixing case. The regions in cyan (pink) indicate where additional 
Higgs bosons close in mass join h = H\ (h = H2) to build up the signal and lead to reduced 
cross sections that differ by more than 10% from the one of h alone. Depending on the value 
of At and the final state these regions are more or less extended: The experimental resolution 
in the various final states is not the same, which has been taken into account by applying a 
different width in the Gaussian smearing of the non h Higgs cross sections, that are added to 
the h final state. Therefore the parameter regions with several Higgs bosons contributing to 
the final state are for WW final states, where the Higgs mass cannot be reconstructed, different 
from the ones for ZZ. The same holds for the fermionic final states. Here the resolution in 
the rr final states is less good than the one in bb, leading to the 'nose' in the plots for h = H\ 
Fig. 11 (left) against the rr final state j^] Another reason for the difference in the extensions 



of the parameter regions is that due to the different Higgs-gauge and Higgs-fermion coupling 



9 The difference in the bb and rr branching ratios due to QCD corrections is small enough not to play a 
significant role here; nor do the negligible At corrections. 
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Figure 10: Reduced cross section into 77 versus the reduced cross section into bb for At = GeV (upper) and 
At — 1 TeV (lower) with h = Hi (left) and h — if 2 (right). Cyan/pink points indicate the signals where at least 
two Higgs bosons with similar masses overlap and the combined reduced cross section deviates by more than 10% 
from the reduced cross section of the individual Higgs boson. Bars: Experimentally measured values with error 
bars (full/ ATLAS, dashed/CMS). 



structures, for a lot of parameter points the non-/i Higgs state contributions to the gauge boson 
final states can not be important enough to induce a change in the rate by more than 10%. This 
is because the Higgs-gauge couplings for small values of tan/3 are dominated by the up-type 
Higgs component. In order to achieve a large enough production for the h Higgs boson through 
gluon fusion its up-type component must be near the SM value, inducing a very small up-type 
component for the other CP-even Higgs bosons due to coupling sum rules, so that they hardly 
decay into massive gauge bosons. The down-type component of the Higgs bosons, however, 
has not been restricted and therefore both the h Higgs boson and the other one(s) with mass 
close by can have equally important couplings to down-type quarks depending on the amount 
of singlet-doublet mixing. 

A substantial amount of scenarios compatible with an excess in the photon final state is hence 
only due to a superposition of Higgs rates stemming from nearly degenerate Higgs bosons. The 
experimental distinction of such scenarios from single Higgs rates, as has been discussed e.g. in 
Ref. [39 1, would be a clear signal of beyond the SM Higgs physics. 

The plots show the strong correlation between the 77 and the massive gauge boson final 
states: In case the increase in the photonic final state is due to an enhanced photon branching 
ratio caused by a suppression in the decay width into bb, this affects the branching ratio into 
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Figure 11: Reduced cross section into 77 versus the reduced cross section into rr for At = GeV (upper) and 
At — 1 TeV (lower) with h = H\ (left) and h — if 2 (right). Cyan/pink points indicate the signals where at least 
two Higgs bosons with similar masses overlap and the combined reduced cross section deviates by more than 10% 
from the reduced cross section of the individual Higgs boson. Bars: Experimentally measured values with error 
bars (full/ ATLAS, dashed/CMS). 



gauge bosons as well and leads also here to larger rates. Should the gauge boson reduced cross 
sections turn out to be exactly SM-like, a strongly enhanced rate into the 77 final state would 
be difficult to comply with. Nevertheless, even in this case enhanced photonic rates up to ~1.6- 
1.8 are still possible. At the present status of experimental errors and experimental resolution 
everything is still compatible. There is a little bit more tension with the CMS results, as CMS 
finds suppressed rates into ZZ, WW contrary to ATLAS reporting enhanced rates. With more 
data accumulated by the experiments and reduced errors on the fixx values future will show 
which of these scenarios will survive and which will be excluded. The correlation between the 
photon and the fermion final states on the other hand is much less pronounced. While in the 
gauge boson final states the branching ratios are simultaneously affected by a change in the bb 
decay mode, the down-type fermion final states are less sensitive to such a change. In the bb 
final state the fj, value reported by ATLAS lies below the allowed regions, the one of CMS above, 
both still compatible within the large errors with the results of the parameter scan so that at 
present no conclusive statement can be made. In the rr final state the reported fi value is below 
one and hence the Higgs-rr coupling suppressed. The ATLAS and CMS values are on the left 
border of the allowed parameter range and compatible within errors, which also in this channel 
are still too large to make firm statements. 
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Concerning perturbativity constraints the regions shown in Figures [8-11 are mostly com- 
patible within the NMSSM without the inclusion of extra matter. A few scenarios require the 
inclusion of extra matter at 1 TeV. In case of tan (3 = 4 the shapes of the parameter regions cor- 
responding to Figures [8 11 stay approximately the same but are much less dense in the amount 
of allowed scenarios. 

With increasing data the precision on the signal strengths reported by the experiments will 
improve and the exclusion limits will become more stringent. This has to be taken into account 
when combining signals stemming from two Higgs bosons which are close in mass. Thereby the 
reduced cross sections fixx wm change. In particular the enhancement in the 77 final state, 
which in a substantial amount of scenarios is due to this superposition, may partially disappear. 
Furthermore, it affects scenarios which have been excluded due to too large signal rates from the 
combination of two Higgs signals. On the other hand, the improved precision would then reveal 
two Higgs signals lying close to each other, providing an unambiguous sign of BSM physics, 
should the NMSSM or some other multi-Higgs sector be realised in nature. 



6 Summary and Outlook 

In this paper we have studied the phenomenology of Higgs bosons close to 126 GeV within the 
scale invariant unconstrained NMSSM, focusing on the regions of parameter space favoured by 
low fine-tuning considerations, namely stop masses of order 400 GeV to 1 TeV and an effective 
\x parameter between 100-200 GeV, with large A (which is required to remain perturbative up 
to the GUT scale) and low tan /3 =2-4. 

By performing scans over the above parameter space, focusing on the observable Higgs cross 
sections into 77, WW, ZZ, bb and rr final states, we have studied the correlations between 
these observables. Although we examined only a limited parameter range in A K , A\ and n, 
we found a substantial amount of parameter space which can lead to Higgs boson masses and 
couplings compatible with the latest LHC results. 

There are basically two types of NMSSM scenarios compatible with the data, corresponding 
to the 126 GeV Higgs boson being either the lightest CP-even Higgs boson Hi or the second 
lightest one Hi- Our results clearly favour the second option, however the first option is still 
possible but it requires additional extra matter at the TeV scale in order to maintain the per- 
turbativity of A, as well as large stop mixing and low tan (3 ~ 2 (for example tan ft = 4 is 
not allowed) and even then the allowed parameter space is relatively sparse. We emphasise that 
these conclusions only apply to the natural NMSSM, in the low fine-tuning region defined above, 
and that larger stop masses and mixing (above one TeV) would allow a larger parameter space 
with Hi at 126 GeV. 

The enhancement in the Higgs rate in the di-photon channel that we observe in our results 
is due to a combination of factors. Firstly there can be an enhancement in the dominant gluon 
fusion Higgs production cross section due to the light squarks in the loop, where light stops 
are a feature of the low fine-tuning region. Secondly the Higgs branching ratio in the di-photon 
channel can be enhanced due to two sub-factors, namely (i) an increased di-photon partial width, 
induced by stop, charged Higgs and chargino loops, and (ii) a suppressed total width, due to 
a suppressed Higgs coupling to b quarks resulting from singlet- doublet mixing. Concerning 
case (i), the stop loop effect in the photonic decay is opposite to the one in gluon fusion and 
depends on the mixing. In the case (ii) also the rates into the fermionic final states can be 
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suppressed. The reported best fit values of the signal strengths [i in the bb and rr final states 
by ATLAS and CMS still suffer from large errors, so that it is difficult to draw a conclusion on 
possibly suppressed couplings to down-type fermions. The allowed parameter ranges we found 
for enhanced di-photon event rates are compatible with the experimental best fits of the fi values 
in the various final states. 

Although our results encompass the SM case where all \i values are equal to unity, we also 
allow for significant and correlated departures from unity for all channels. While we do not 
find any significant correlation between di-photon and fermion rates (fx values), we do find a 
correlation between the di-photon and massive gauge boson rates (// values). However, given 
the present status of the experimental accuracy in the various final states, it is not possible 
to draw any conclusions about this. Nevertheless it is clear that the natural NMSSM Higgs 
sector (corresponding to the low fine-tuning region as defined above) is nicely compatible with 
all experimental results, with the bulk of the data points corresponding to the the second lightest 
Higgs boson H2 having a mass of about 126 GeV. In this favoured case, the 126 GeV H2 boson 
can decay into pairs of lighter neutralinos, CP-even or CP-odd Higgs bosons, providing a smoking 
gun signature of the NMSSM. 

We also emphasise that a good part of the parameter space involves a Higgs spectrum where 
the two lighter CP-even Higgs bosons H\ and H2 are close in mass. It is then the combination 
of their reduced cross sections and rates which is observed in the experiment. However, with 
increasing accuracy in the Higgs boson mass resolution, future LHC data may resolve these two 
states. Observing the two separate CP-even Higgs bosons Hi and H2 with different masses 
would not only rule out the Standard Model, but could also provide direct evidence for the 
Higgs sector of the NMSSM. 

With future LHC results, the best fit values of the signal strengths and their errors in the 
different channels will change, leading to different positions and error bars on the data points 
represented by crosses in our plots. However, the overall pattern of the plots themselves will 
not change substantially. Thus future data can be compared to our predictions to check the 
compatibility of the natural NMSSM with experiment. If stops are not discovered below one 
TeV, and instead the experimental limit on the stop masses increases, then the range of the stop 
masses and mixing may need to be extended beyond the low fine-tuned region considered here, 
leading to enlarged parameter regions of the NMSSM. However, for the moment, the natural 
NMSSM is still viable, with the characteristic Higgs spectrum and properties discussed in this 
paper. 
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Appendix 



A Best fit values of the signal strength 

We list the best fit values of the signal strengths fi = cr/asM in the various final states reported 
by ATLAS 



64 



66 



and CMS (2 67 -70 , which we have applied in our plotsQ 



Experiment 



ATLAS 



CMS 



Final state 



s, L) 



77 

WW 

ZZ 

bb 

TT 



77 
WW 

zz 
bb 

TT 



(7 TeV, 4.8 fb- 1 )+(8 TeV, 5.9 ftr 1 ) 

(8 TeV, 13 fb- 1 ) 
(7 TeV, 4.8 fb- 1 )+(8 TeV, 5.8 fb -1 ; 
(7 TeV, 4.7 fb -1 )+(8 TeV, 13 fb" 1 ) 
(7 TeV, 4.6 fb _1 )+(8 TeV, 13 fb^ 1 ) 



(7 TeV, 5.1 fb- 1 )+(8 TeV, 5.3 fb" 1 ) 

(7 TeV, 4.9 fb- 1 )+(8 TeV, 12.1 fb" 1 ; 
(7 TeV, 5.1 fb- 1 )+(8 TeV, 12.2 ft) -1 ) 
(7 TeV, 5 fb" 1 )+(8 TeV, 12 fb" 1 ) 
(7 TeV+8 TeV, 17 fb" 1 ) 



\i = a/asM 



64 



1.8 ±0.5 1 
1.5 ±0.6 
1.4 ±0.6 |] 
-0.4 ± 1.1 
0.7 ±0.7 66 



65 



i.56 ± (uair 

0.74 ± 0.25 |67| 

,+0.35 



0.8 
1.3 



0.28 

-0.7 
-0.6 



68 



69 



0.72 ± 0.52 1701 



Table 1: Best fit values of the signal strength fi in the 77, WW, ZZ, bb and rr final states reported by ATLAS 



1 64 66 and CMS 2 67 70 



10 The statistical (±0.7) and systematic error (±0.8) in the bb final state reported by ATLAS have been added 
in quadrature. 
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